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From  detailed  analyses  of  both  data  sets,  we  conclude  that  IPS  prediction 
of  the  onset  of  geomagnetic  disturbances  may  be  feasible  with  a lead  time 
of  about  one  day.  . 


ru:ssi3«l»  1 

1 »7it 

wmtf  SKtiM 

I«t1  SKllM 

□ 



: tTKisniai/usiuBiui’t  caoa 
01:1  hUL  ui~ir tfeCIAl 


D D C 

rrvrT-r^'^-n  ns 

U<i  FEB  21  1978 

JljEisEnirE 

D 


Unclass i f led 


Classification  of  this  FAc.E'*Fh»f»  oaia 


1 


INTRODUCTION 


We  proposed  a 15-month  Investigation  of  the  feasibility  of  using 
Interplanetary  radio  scintillation  (IPS)  observations  at  34.3  MHz  for  the 
prediction  of  solar  particle  events  and  geomagnetic  activity.  Our  overaLt 
oonalusion  is  that  24-hour  IPS  prediction  of  eolaop-terreatrial  diaturbanoee 
may  be  feaaible  within  certain  conatrainta  intrinaia  to  the  IPS  technique. 
Since  the  observations  upon  which  we  base  our  conclusions  were  limited  to 
two  periods,  May-December  1974  and  May-August  1976,  the  Inferences  drawn 
directly  from  our  data  are  strictly  applicable  only  to  the  descending  phase 
of  the  solar  cycle.  However,  our  theoretical  analysis  allows  us  to  make 
some  evaluation  of  IPS  prediction  feasibility  during  the  rise  and  may 1 mum 
of  the  solar  cycle. 

A promising  new  development  Is  evidence  that  Important  additional 
prediction  Information  Is  contained  In  the  functional  form  of  the  IPS  power 
apeatrtn.  Previous  comparisons  of  IPS  observations  with  interplanetary 
structure  relied  mainly  on  the  IPS  Index  which  Is  essentially  the  total 
normalized  power  Integrated  over  the  IPS  po%ier  spectrum. 

This  report  will  not  cover  the  actual  operations  of  the  COCOA-Cross 
array  by  the  University  of  Iowa  at  the  Clark  Lake  Radio  Observatory.  They 
were  sumnarlzed  in  our  Interim  Scientific  Report  (1  April  1976-31  March  1977) 
on  Contract  NPP-75-157  from  the  Air  Force  Office  of  Scientific  Research  who 
are  also  supporting  our  IPS  prediction  effort.  Progress  on  the  up-grading 
of  the  array,  made  necessary  by  the  flooding  of  Clark  Dry  Lake  by  Hurricane 
Kathleen  In  the  Fall  of  1976,  was  again  delayed  by  an  unprecedented  15-lnch 
flooding  In  August  1977.  Consequently,  the  bulk  of  our  effort  under  this 
contract  was  devoted  to  the  scientific  analysis  IPS  synoptic  Index  data 
taken  May-December  1974  by  the  University  of  Iowa,  and  to  the  IPS  synoptic 
power  spectra  taken  May-August  1976  with  the  digital  data  acquisition  and 
reduction  systems  developed  at  the  Applied  Physics  Laboratory.  The  1976 
OOCOA-Cross  data  were  significantly  supplemented  by  concurrent  observations 
using  the  University  of  Maryland  TPT  array  at  CLRO,  for  which  we  are  Indebted 
to  Professor  W.  C.  Erickson. 

The  results  of  the  scientific  analysis  provided  considerable  Insight 
Into  the  prediction  problem  and  were  submitted  as  three  separate  articles  to 


the  Journal  of  Geophyaiaal  Reaearch.  All  three  have  been  refereed  and  the 
revisions  are  pending  acceptance  by  the  editor.  Many  of  the  Ideas  sunnnar- 
ized  In  this  Report  are  discussed  more  fully  In  these  papers  which  are 
currently  available  in  their  preprint  form. 

Interplanetary  Sointillation  at  Large  Elongation  Anglea:  Reaponae 
to  Solar  Vind  Denaity  Structure,  by  F.  T.  Ersklne,  W.  M.  Cronyn,  S.  D.  Shawhan, 
E.  C.  Roelof  and  B.  L.  Gotwols  (University  of  Iowa  Preprint  77-24). 

Analyaia  of  Interplanetary  Sointillation  Spectra  at  Large  Elongation 
Anglea,  by  D.  G.  Mitchell  (Applied  Physics  Laboratory/ JHU  Preprint  77-07). 

Synoptic  Analyaia  of  Interplanetary  Radio  Scintillation  Spectra 
Obaerved  at  34  MHz,  by  B.  L.  Gotwols,  D.  G.  Mitchell,  E.  C.  Roelof,  W.  M. 
Cronyn,  S.  D.  Shawhan  and  W.  C.  Erickson  (Applied  Physics  Laboratory/ JHU 
Preprint  77-13) . 

2.  THEORETICAL  BACKGROUND 

One  of  the  Justifications  for  constructing  the  COCOA-Cross  array  to 
operate  at  34.3  Mlz,  a frequency  lower  than  usually  utilized  for  IPS  observa- 
tions, was  that  the  increased  IPS  response  at  the  lower  frequency  would  make 
it  possible  to  observe  sources  at  large  elongation  angles  (c)  from  the  sun 
(where  IPS  response  falls  off  due  to  both  the  decreased  mean  density  of  inter- 
planetary electrons  and  the  reduction  in  the  solar  wind  velocity  component 
transverse  to  the  line  of  sight  to  the  source).  The  expectation  was  to  have 
a dense  grid  of  scintillating  sources  covering  the  sky  above  the  radio  horizon 
(except  for  the  region  occupied  by  the  galactic  plane  where  interstellar  scat- 
tering reduces  the  response  to  interplanetary  scattering). 

The  main  disadvantage  in  going  to  a lower  observing  frequency  lies 
in  the  increase  in  the  apparent  angular  diameter  of  the  source,  due  to  the 
increase  with  wave  length  of  both  interstellar  scattering  and  extent  of 
extended  sources.  The  apparent  angular  source  diameter  is  an  important  para- 
meter in  IPS  prediction,  since  it  determines  from  how  far  out  into  the  inter- 
planetary medium  detectable  scintillation  power  is  generated.  To  be  more 
•xact,  in  weak-scattering  theory  the  frequency  (v)  dependence  of  the  contri- 
bution to  the  IPS  power  spectrum  dP(v)  from  an  interval  dL  along  the  line  of 
sight  contains  a source  visibility  function  which  we  may  take  in  the  form 

exp(-v^/vc^) , where  v_  - V /II0L.  Here  0 is  the  apparent  angular  diameter  of 
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the  source  and  is  the  component  of  the  solar  wind  velocity  transverse  to 
the  line  of  sight.  For  a given  frequency  v In  the  IPS  power  spectrum,  contri- 
butions from  distances  L > V^/TTev  will  be  strongly  attenuated.  We  measure 
P(v)  over  a range  0.1  < v < 1.5  Hz,  and  apparent  angular  diameters  are  often 
> 1 arc-second  at  34.3  MHz.  Using  a typical  value  “■  400  km/s  and  v “ 0.5 
Hz,  the  nominal  region  "probed"  by  this  component  of  the  IPS  spectrum  for  a 
1"  source  Is  within  0.35  AU  of  the  earth.  Of  course.  If  the  Interplanetary 
turbulence  were  localized  beyond  0.35  AU,  It  would  still  contribute  at  v * 0.5 
Hz  If  It  were  strong  enough. 

There  Is  another  Important  process  Intrinsic  to  IPS  which  affects  the 
power  spectrum.  This  Is  "Fresnel  filtering"  and  It  Is  Independent  of  source 
size.  The  scintillation  power  contributed  by  electron  density  Irregularities 
at  a distance  L rises  with  decreasing  scintillation  frequency  v,  but  becomes 
constant  below  the  "Fresnel  frequency"  ~ This  occurs  because 

Irregularities  whose  scale  sizes  exceed  the  diameter  of  the  first  Fresnel 
Interference  zone  are 'Ylltered"  out  of  the  spectrum  by  destructive  Inter- 

ference. Thus  Irregularities  near  the  earth  are  less  effective  at  contribu- 
ting power  at  low  scintillation  frequencies  than  those  further  away.  For  a 
scintillation  frequency  v,  Fresnel  filtering  will  occur  for  all  distances 
L < /HXv  . At  X ■ 10  m observing  wavelength  and  “ 400  km/s,  this  gives 
Fresnel  filtering  for  L <0.14  AU  for  v ■ 0.5  Hz.  Note  that  as  one  goes  to 
lower  observing  frequencies  (larger  X),  the  Fresnel  effect  Is  reduced. 

These  two  effects  were  combined  In  a theoretical  calculation  of 
the  thin  screen  power  spectrum  by  iHtchell  and  Roelof  (J.  Geophye.  Rea., 

81 j 5071,  1976).  The  contribution  at  each  frequency  was  then  Integrated  using 
all  screens  along  the  line  of  sight  for  a spherically  symmetric  medium.  The 
spectra  are  shown  In  Figure  1 and  demonstrate  the  Fresnel  flattening  at  low 
scintillation  frequencies  and  the  source  size  suppression  at  higher  frequencies. 

It  Is  clear  therefore  that  the  form  of  the  observed  IPS  power  spectrum 
contains  Information  on  the  distribution  of  turbulence  within  the  medium,  since 
these  two  filtering  effects  on  the  power  spectrum  must  be  weighted  by  the  varia- 
tion for  the  distribution  of  the  turbulence  along  the  line  of  sight.  Even  In  a 
quiescent  Interplanetary  medlua,  the  turbulence  distribution  Is  a strong  func- 
tion of  the  elongation  angle  c from  the  sun  to  the  source.  Figure  2,  taken  from 
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Mitchell  (1977)  shows  the  theoretical  contribution  dm/dL  along  the  line  of 

sight  for  the  band-passed  power  spectrum  m (0.1-1. 5 Hz)  observed  at  34.3  MHz 

from  a spherically  symmetric  Interplanetary  medium  with  the  electron  variance 
2 -4 

< 6N  > varying  as  r with  the  distance  r from  the  sun.  The  curves  cover  the 
range  of  elongations  60**  < e < 160*  and  were  calculated  for  source  ant.ular 
sizes  e ■ 0.01"  (effectively  a point  source),  6 ■ 3.0"  (not  an  uncommon  size  at 
low  frequencies)  and  6 ••  7.0"  (a  large  source).  The  Fresnel  filtering  Is 
apparent  for  small  distances  L,  while  the  source  size  filtering  restricts  the 
quiescent  medium  response  to  ^ 0.5  AU.  At  34.3  MHz,  scintillation  Is  depressed 
or  "quenched"  at  elongations  e < 40“  because  the  scattering  is  strong  due  to 
the  high  electron  density. 

2 

Although  these  calculations  assumed  a quiescent  medium  with  < 6N  > 

-4 

r , they  can  be  generalized  to  a disturbed  medium  with  locally  enhanced  tur- 
bulence. Mitchell's  calculations  for  Figure  2 actually  divided  the  medium  Into 
20  slabs,  and  the  locations  of  slabs  7 and  13  (Indicated  on  Figure  2)  are  plotted 
along  lines  of  sight  at  e ■ 60“  and  120“  In  the  ecliptic  plane  In  Figure  3,  along 
with  a sketch  of  a corotating  turbulent  region  (shovm  fixed  for  sequential  posi- 
tion of  the  earth  every  two  days,  a-d). 

3.  IMPLICATIONS  OF  THEORY  FOR  PREDICTION 

Ue  now  have  the  essential  concepts  to  Interpret  the  COCQA-Cross  obser- 
vations. Referring  back  to  Figure  2,  we  see  that  the  bulk  of  the  response  of 
the  3"  source  Is  within  slabs  1-7.  Therefore  for  the  geometry  of  the  aorotating 

turbulent  region  shown  In  Figure  3,  enhanced  scintillation  at  e **  60“  would  be 

expected  not  more  than  1-2  days  before  the  corotating  turbulence  reached  earth 
(c),  and  would  cease  about  a day  afterwards  (b).  No  IPS  enhancement  would  be 
seen  at  c ■ 120“  much  In  advance  of  the  arrival  of  the  turbulence  (c) , but  the 
enhancement  may  last  about  a day  longer  (b)  than  that  at  60“  (unless  the  tur- 
bulence weakens  outside  the  earth's  orbit,  as  It  well  may).  Most  interplanet- 
ary disturbances  during  1974-6  were  recurrent  and  corotating,  and  we  shall  show 

below  that  these  anticipated  patterns  occurred  In  the  COCOA-Cross  data.  For 
radially  propagating  flare-associated  dlsttirbances,  the  response  pattern  would 
be  similar,  but  of  somewhat  shorter  duration,  as  can  be  seen  by  sketching  an 
expanding  spherical  shell  centered  on  the  sun  on  Figure  3. 


The  theory  tells  us  that  the  enhancement  of  total  Bcintillation 

porcer  is  not  the  sole  predictor  available  for  the  arrival  of  turbulence. 

Referring  to  Figure  1,  we  see  that  the  shape  of  the  power  speatrum  itself 

will  change  as  the  turbulent  region  approaches.  For  the  -g"  source  at  e = 

60°  in  Figure  lb,  it  may  be  seen  that  the  contribution  from  the  more  distant 

regions  is  narrower  (due  to  the  finite  source  size  filtering  v > v “ L ^) , 

S 

whereas  the  contributions  from  the  closer  region  is  broader  (because  the 
source  size  filtering  is  less  important  but  the  Fresnel  filtering  flattens 
the  spectrum  for  v < <x  L ^) . We  have  found  this  spectral  broadening 
diagnostic  in  the  data,  and  its  Importance  lies  in  the  additional  informa- 
tion on  distance  to  the  turbulent  region  which  it  contains. 

4.  SYNOPTIC  OBSERVATIONS;  SCINTILLATION  INDEX 

Data  Analysis 

The  normalized  scintillation  index  curves  for  a cluster  of  six 
sources  covering  an  ecliptic  latitude  range  of  9°  and  a longitude  range  of 
57°  are  displayed  in  Figure  4.  Over  the  time  period  1974,  days  231  to  296, 
this  cluster  was  west  of  the  sun.  IMP-7  and  8 (near  the  Earth)  solar  wind 
density  and  velocity  measurements  (from  the  Los  Alamos  Scientific  Laboratory 
detectors),  the  planetary  Ap  index  and  the  inferred  interplanetary  magnetic 
field  polarity  at  the  Earth  are  also  provided.  There  is  good  temporal 
agreement  between  the  occurrence  of  density  enhancements  (which  are  often, 
but  not  always,  associated  with  the  magnetic  sector  boundaries  and  which 
are  followed  by  velocity  peaks),  and  the  occurrence  of  Increases  in  the 
scintillation  index  for  each  of  the  sources  individually.  This  correlation 
is  expected  because  of  the  Increase  in  scattering  associated  with  the  density 
increase  along  the  line  of  sight.  The  density -as  so  dated  IPS  index  enhance- 
ments precede  large  increases  in  geomagnetic  activity  (Ap)  by  1-2  days.  How- 
ever, the  quantitative  correlation  of  index  response  from  source-to-source 
varies.  For  the  enhancement  on  days  256-259  all  sources  exhibit  approximately 
the  same  Increase  in  index.  For  days  231-234,  however  3C68.2  barely  responds 
in  comparison  to  the  other  sources.  Also,  3C125  and  3C153  Increase  signifi- 
cantly less  on  days  285-289  than  the  other  sources.  These  observations  of 
distinct  source-to-source  differences  within  a source  cluster  suggest  spatial 
variations  in  the  electron  density  turbulence  on  the  scale  of  2-10°  and/or 
temporal  variations  on  a scale  of  3-16  hours.  This  scale  size  for  density 
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structure  is  similar  to  velocity  structure  deduced  for  the  solar  wind  velocity 
difference  between  two  spatially  separated  satellites  where  a longitudinal 
coherence  time  of  18  hours  (10®)  is  found,  as  well  as  significant  changes  on 
the  scale  of  a few  degrees  in  latitude. 

It  is  apparent  from  Figure  4 that  missing  observations  obscure  the 
pattern  over  the  sky  on  a given  day.  A further  difficulty  is  encountered  when 
we  desire  to  Intercompare  different  sources  on  a given  day;  how  are  the  IPS 
responses  of  different  sources  to  be  normalized?  In  Figure  4 the  normalization 
to  peak  values  is  really  arbitrary,  since  there  is  no  assurance  that  each  source 
line-of-sight  contained  the  same  turbulence  on  the  peak  days.  Therefore  we 
developed  a non-parame tr ic  ranking  scheme  and  also  grouped  adjacent  sources 
in  clusters  (to  minimize  the  effects  of  missing  observations).  We  selected 
45  of  the  most  reliable  scintillating  sources  for  a source  "clustering  analysis". 
As  shown  in  Figure  5,  we  chose  10  regions  or  "boxes"  in  the  sky,  each  contain- 
ing up  to  8 radio  sources,  with  the  sources  spaced  no  more  than  about  30®  from 
any  other  source  in  the  same  box.  The  10  sky  boxes  span  nearly  all  ecliptic 
longitudes.  Most  of  the  boxes  are  centered  at  mid-northern  ecliptic  latitudes 
(from  20®  to  50®).  However,  Box  5 is  just  south  of  the  ecliptic,  while  Box  3 
(the  single  source  3C161)  is  at  -29®  latitude.  Box  9 is  centered  at  a rather 
high  latitude  of  about  70®.  The  method  of  assigning  a weighted  IPS  response 
to  each  box  (Qj^)  is  described  in  detail  by  Erskine  et  al.  (1977). 

With  the  set  of  weighted  box  responses  Qj  ( i = 1,  ...,  10)  for  each 
day  of  observations,  we  can  now  study  the  general  large  scale  pattern  of  scin- 
tillation over  the  sky.  One  graphical  representation  of  the  is  given  for 
the  three  observational  periods  in  Figures  6,  7 and  8.  Four  levels  of  shading 
are  used  in  Figures  6-8  to  display  the  high,  moderate,  low  or  missing  IPS 
activity  estimates  for  each  sky  box,  each  day,  with  the  darkest  shading  indi- 
cating "high",  and  no  shading  indicating  "missing  observation".  The  boxes 
are  arranged  from  top  to  bottom  in  order  of  decreasing  ecliptic  longitude. 

At  the  bottoms  of  Figures  6-8  are  displayed  correlative  interplanetary  and 
geomagnetic  data  in  a format  similar  to  Figure  4 (see  captions  for  Figures 
6-8).  The  solar  wind  measurements  are  from  the  Los  Alamos  Scientific  Labora- 
tory (lASL)  detector  on  IMP  7 and  8,  and  the  magnetic  field  from  the  Ihilverslty 
of  Iowa  plasma  wave  search  coil  detector  on  IMP-8  spacecraft.  In  addition,  we 
have  computed  the  average  "all-sky  IPS  acticity  index",  (Q) , over  all  boxes. 
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We  also  compute  ^^^west  ®3ch  day,  using  only  the  eastern  and 

western  sky  boxes.  We  show  <Q>,  <^>west  ^‘^^high  lat.  ‘^9 

time  at  the  centers  of  Figures  6-8. 

Examination  of  Figures  6-8  for  evidence  of  a spatial-temporal  pattern 
of  IPS  activity  reveals  the  presence  of  vertical  bands  of  occasional  very  low 
IPS  acticity  and  "events"  of  high  and  moderate  IPS  activity.  These  event  bands 
encompass  nearly  all  ecliptic  longitudes  and  appear  to  be  most  prominent  near 
times  of  solar  wind  density  Increases,  The  plots  of  ^Q),  ^^^v,est 

all  peak  at  approximately  the  same  time.  The  peaks  of  IPS  activity 
appear  to  coincide  well  with  peaks  in  solar  wind  proton  density,  and  for  the  most 
part,  1 B 1 also.  The  peaks  in  geomagnetic  Ap  occur  sei'eral  days  later  than 
the  peaks  in  IPS  activity. 

Correlation  Analysis 

In  order  to  examine  the  temporal  relationship  between  the  34.3  MHz 
IPS  enhancements  and  correlative  ground-based  geomagnetic,  solar  disc  observa- 
tions, and  in  situ  spacecraft  solar  wind  data,  Erkslne,  in  his  dissertation, 
performed  a time  series  cross-correlation  analysis  between  < Q(t)  > and  the 
correlative  variables.  The  data  base  for  this  cross-correlation  spans  1974, 
t »=  day  150  to  t = day  362.  The  range  of  time  shift  used  is  -10  days  s x s 
159  days,  corresponding  to  75%  of  the  data  base.  < Q(t)  > is  the  all-sky  IPS 
activity  index  as  described  earlier  and  Q is  the  mean  < Q(t)>  over  the  data 
base.  The  correlation  has  been  performed  against  the  following  parameters 
X(t):  all-sky  IPS  activity  measure  Itself  (auto-correlation)  < Q >,  < Q > , 

< Q ^ Q lat  * wind  proton  density  N,  interplanetary  magnetic 

field  magnitude  (ecliptic  component)  | B | , solar  wind  proton  bulk  velocity 

xy 

V,  geomagnetic  planetary  A index,  occurrences  of  Interplanetary  sector  crossings, 
geomagnetic  storm  sudden  commencements.  Type  II  and  IV  radio  bursts  (on  the 
eastern  and  western  solar  hemispheres),  as  well  as  central  meridian  passage  of 
coronal  "hole"  western  edges  as  estimated  from  Figure  3 of  Hansen  et  al.  (Planet. 
Space  Soi.,  24,  381,  1976).  In  order  to  compute  a cross-correlation  function 

for  the  discrete  events  (SSC's  Type  II  and  IV  bursts  and  coronal  hole  CMP's), 
each  event  was  assigned  a triangular  time  profile  centered  on  that  date  with 
a width  of  one  day  at  half-maximum. 
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Six  representative  correlation  functions  are  presented  in  Figure  9. 
Above  each  plot  of  r(T)  we  have  indicated  the  time  shifts  (days)  for  most  of 
the  peaks.  The  value  of  maximum  correlation  coefficient  nearest  t ^ 0 and  the 
number  of  points,  N,  in  the  time  series  are  also  indicated.  The  Inclined  lines 
are  99Z  confidence  limit  estimates  (assuming  each  measurement  is  an  Independent 
estimate) . Missing  observations  due  to  radio  frequency  Interference  and  tele- 
scope "down  time"  account  for  the  fact  that  N is  considerably  less  than  the 
number  of  days  spanned  by  the  data  base.  We  also  summarize  in  Table  1 the 
essential  parameters  of  all  thirteen  cross-correlation  functions  which  we  com- 
puted: the  amplitudes  and  lags  of  the  correlation  peaks  nearest  to  zero  lag, 

and  "average"  lag  which  is  the  mean  lag  (modulo  27  days)  over  all  the  identi- 
fiable recurrence  peaks. 

We  note  that  99%  confidence  limit  estimates  of  r = 0.21  fo’^  n = 135 
assume  that  each  day's  observation  constitutes  an  Independent  measurement.  This 
is  not  true  for  IPS,  solar  wind  and  geomagnetic  disturbances  in  which  "events" 
last  3 days.  This  is  the  well-known  conservation  effect  in  the  statistical 
analysis  of  geomagnetic  data.  The  net  result  is  that  the  number  of  "independent" 
points  will  then  be  approximately  one-third  of  the  number  of  observations,  so 
that  the  99%  limit  for  n =“  135  (although  applicable  to  discrete  events  such  as 
sector  boundaries,  SSC's,  radio  bursts,  and  coronal  hole  edge  CMP's),  should  be 
interpreted  as  a 90%  confidence  limit  with  n 45  for  the  persistent  enhance- 
ments of  IPS,  solar  wind  velocity  and  density,  IMF  intensity  and  Ap  index. 

We  see  from  Figure  9 that  the  IPS  enhancements  recur  consistently 
over  the  Interval  spanned  by  the  time  shifts  (5  solar  rotations) . The  autocor- 
relation of  < Q > has  a mean  recurrence  period  of  28  ± 2 days,  verifying  our 
interpretation  that  the  dominant  IPS  structure  during  this  period  is  corotating. 

The  cross-correlations  of  < Q > against  <Q>  ^ and  < Q ^ 

east  west  high  lat. 

(not  shown)  give  mean  recurrence  periods  of  27  ± 2 days,  28  ± 2 days  and  27  ^ 1 
day  respectively.  There  is  no  significant  time  shift  between  the  correlations 
with  eastern,  western  and  high  latitude  sources  (to  within  an  uncertainty  of 
1-2  days),  and  the  correlation  coefficients  in  Table  1 are  quite  high  ( >0.76 
for  T " 0 days).  Many  of  the  correlations  in  Figure  9 show  subsidiary  peaks 
located  11-15  days  from  the  main  peaks  which  are  approximately  at  t ■ k • 27 
days,  k ■ 0,  1,  2,  ...,  5,  although  they  usually  fall  below  the  99%  confidence 
limit  for  k a 2.  The  subsidiary  peaks  (which  are  occasionally  bifurcated)  in 
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the  correlation  plots  of  Figure  9 are  the  result  of  the  somewhat  asymmetric 
solar  wind  two-sector  structure  in  1974.  Computer  simulations  have  shown 
that  the  amplitude  of  the  subsidiary  peak  in  r(T)  is  far  more  sensitive  to 
changes  of  spacing  of  the  two-sector  structure  than  to  changes  in  height  of 
the  maxima. 

Further  examination  of  Figure  9 and  Table  1 confirms  our  impression 
from  Figures  6-8  that  the  IPS  actlcity  index  < Q > correlates  better  with  solar 
wind  proton  density,  r(0)  •=  0.56,  than  with  any  of  the  other  correlative  para- 
meters examined,  and  that  density  and  IPS  activity  are  evidently  in  phase  with 
each  other.  IPS  activity  is  also  in  phase  with  interplanetary  magnetic  field 
magnitude,  but  the  correlation  is  noticeably  smaller  r(0)  = 0.37.  As  can  be 
seen  from  Figures  6-8,  most  large  increases  in  | B 1 are  the  result  of  density 
compressions  in  stream-stream  interactions.  The  IPS  activity  anticorrelates 
well  with  solar  wind  proton  bulk  velocity,  r(-2)  = -0.52,  as  well  as  correlating 
at  a later  lag  r(+3)  = 0.38.  These  are  the  normal  lags  between  velocity  minima 
preceding  density  peaks  and  velocity  maxima  following  them,  as  is  apparent  in 
Figures  6-8. 

The  correlations  with  geomagnetic  activity  are  also  consistent  with 
the  nominal  structure  of  stream-stream  interactions.  The  correlation  with  SSC's 
peaks  at  zero  lag  while  peaks  a day  later  (as  the  storm  develops),  but  with 
a correlation  coefficient  r(+l)  • 0.43  that  is  higher  than  for  the  discrete  SSC's. 

The  weakest  correlations  obtained  were  between  < Q > and  IMF  sector 
crossings,  radio  bursts,  and  coronal  hole  CMP’s.  In  appearance,  the  correlation 
functions  are  rather  similar  and  are  marginally  statistically  significant,  so 
only  the  result  for  sector  crossings  appears  in  Figure  9.  The  sector  crossings 
would  be  expected  to  peak  at  zero  lag,  due  to  their  statistical  associations 
with  SSC's,  but  the  correlation  peaks  for  Type  II  and  IV  bursts  and  coronal 
holes  are  not  well  enough  defined  to  justify  further  interpretation,  although 
in  our  judgement  the  dominant  class  of  IPS  disturbances  during  May-December 
1974  exhibited  a high  degree  of  recurrence  and  were  not  flare-associated. 

Event  Analysis 

We  pointed  out  above  that  the  density  events  are  sometimes  preceded  by 
a ttough  of  moderate  and  low  IPS  activity.  It  seems  that  there  was  a compara- 
tively non-turbulent  region  west  of  the  density  structure  which  allowed  it  to 
be  defined  more  clearly  by  line-of-sight  measurements.  Thr  trough  most  evident 
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in  Figure  6 is  that  in  the  eastern  sources  (boxes  5-8)  on  days  156-159.  Even 
though  we  are  limited  by  Incomplete  coverage  and  somewhat  noisier  data  during 
this  Initial  period  of  our  observations,  the  data  still  permit  a more  detailed 
examination  of  the  cluster  analysis  for  the  duration  of  the  event,  days  156-163. 

Plotted  In  the  left-hand  column  of  Figure  10  are  the  for  each  day 
when  observations  were  available.  Since  3C161  was  not  observed  at  this  time, 
there  are  no  data  for  Q^.  At  the  head  and  foot  of  the  column  are  the  LASL  solar 
wind  densities  (3-hour  averages).  The  orientation  of  the  sky-boxes  relative  to 
the  sun  Is  shown  In  the  right  half  of  the  figure.  Boxes  5-7  display  low  Q-values 
on  days  156-158  while  Box  4 Is  low  on  day  157  but  shows  a moderate  enhancement 
(Q  = 0.9)  on  day  158.  Box  4 would  be  the  first  group  of  sources  to  detect  a 
corotating  structure.  The  eastern  onset  of  the  event  Is  clear  on  day  159  In 
Boxes  5 and  6 (no  observations  for  Box  4) , while  Boxes  7 and  8 (covering  solar 
ecliptic  longitude  between  90°  and  150°),  remain  at  Q ■ 0.  The  solar  wind 

_3 

density  peaks  at  N - 25  cm  early  on  day  161,  and  an  "all-sky"  event  occurs 
on  day  162.  The  enhanced  Q- Indices  are  sustained  through  day  163  in  Boxes  7-10 
(viewing  sources  at  longitude  west  of  the  direction  of  a stream  front  at  earth) , 
while  Box  1,  the  closest  western  box  to  the  sun,  shows  an  Increase  on  that  day. 

On  the  contrary.  Boxes  2 and  4,  just  east  of  the  sun,  have  dropped  to  Q ■ 0 on 

day  163,  suggesting  that  the  disturbance  now  appearing  In  the  west  has  passed 

over  the  earth  and  Is  no  longer  visible  east  of  the  sun.  A sketch  of  a density 

structure  consistent  with  these  observations  Is  offered  In  Figure  10. 

There  Is  an  Interesting  enhancement  In  the  western  Box  1 on  days  156 
and  157.  At  that  time  McMath  plage  region  12972  was  In  the  western  hemisphere 
on  the  sun,  producing  an  Importance  IB  flare  on  day  155.0  (S14,  W13)  and  a 3B 
flare  on  day  157.7  (S16,  W48) . Both  flares  had  accompanying  Type  II  and  IV 
radio  bursts.  It  may  well  be  that  the  sources  In  Box  1 (which  Is  at  low  lati- 
tude) were  responding  to  plasma  ejected  from  the  vicinity  of  MPR  12976  (as 
Indicated  on  the  accompanying  sketch  In  Figure  10) . 

We  shall  now  turn  to  another  IPS  event  wherein  almost  Identical  cir- 
cumstances are  found.  However,  since  this  second  event  occurred  half  a year 
later,  our  source  distribution  favored  detailed  analysis  of  structure  to  the 
west  of  Earth,  with  minimal  coverage  to  the  east.  Figure  11  presents.  In  an 
Identical  format  to  Figure  10,  the  IPS  density  event  of  days  330-338.  Source 
clusters  In  Boxes  6-8  now  view  sources  to  the  west  of  the  sun,  and  they  all 
detect  an  IPS  enhancement  on  days  332  and  333.  At  this  time,  McMath  plage 
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region  13343  was  west  of  central  meridian,  later  producing  on  day  334.3  an 
Importance  -B  flare  (N03,  W52)  accompanied  by  an  Isolated  Type  II  burst.  The 
confirmation  by  three  source  clusters  of  the  western  event  leads  us  to  Inter- 
pret this  western  IPS  enhancement  as  plasma  emitted  from  MPR  13343  during  Its 
activity  prior  to  the  flare  on  day  334.3. 

If  we  now  turn  our  actentlon  to  the  sole  source  cluster  viewing  east 
at  mid-latitudes  (sky-box  10),  there  Is  a clear  enhancement  on  days  333  and  334 
following  low  scintillation  activity  on  days  330-332  and  leading  to  sustained 
high  scintillation  on  days  335-337  which  drops  abruptly  on  day  338.  This 
sequence  of  activity  corresponds  nicely  to  the  solar  wind  density  structure 
that  swept  over  the  earth  on  days  335-337.  It  presumably  would  have  been 
detected  approaching  from  the  east  on  days  333  and  334  and  could  no  longer  be 
viewed  to  the  east  on  day  338  after  passing  the  earth. 

The  western  and  eastern  signatures  are  made  even  more  significant 
by  the  extremely  quiet  conditions  In  the  other  source  clusters  (Boxes  1-5) 
throughout  days  330-334.  It  Is  not  until  days  336  and  337  that  these  sources, 
all  at  elongation  angles  exceeding  90",  fully  respond  to  the  solar  wind  density 
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structure  which  peaked  at  N ••  51  cm  late  on  day  336.  The  suggested  interplan- 
etary density  structures  are  sketched  in  Figure  11.  It  seems  clear  that  It  was 
only  because  the  Interplanetary  medium  to  the  west  of  the  corotating  structure 
had  been  so  quiet  (examine  the  "trough"  which  stands  out  so  clearly  In  Figure  8 
for  days  330-334),  that  the  signatures  of  the  corotating  event  In  the  east  and 
the  active  region  event  In  the  west  were  so  easily  separable.  One  therefore 
Is  led  to  regard  the  less  ordered  structures  seen  In  Figures  6-8  as  perhaps 
reflecting  more  accurately  the  usual  complexity  of  the  density  structure  of  the 
Interplanetary  medium. 

As  a final  detailed  example  of  IPS  events  viewed  over  a wide  range 
of  elongation  angles,  let  us  look  at  the  resporte  of  each  source  observed  on 
day  332  (the  western  event)  and  days  336  and  337  (the  all-sky  density  event) . 

We  present  In  Figure  12  an  equal-area  projection  of  the  sky  as  viewed  from 
the  earth  looking  toward  the  sun.  The  IPS  response  q^  for  each  source  appears  as 
an  open  circle  for  low  activity  (0  £ q^  < 2/3),  a cross  for  moderate  activity 
(2/3  & qj  & 4/3)  and  a star  for  high  activity  (4/3  < q^  £ 2) . The  active 
region  event  Is  evident  on  day  332  (even  though  It  Involves  only  the  moderate 
enhancements  expected  at  elongation  angles  < 60"  due  to  quenching  by  strong 
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scattering).  It  is  obvious  that  this  pattern  of  activity  would  have  gone 

undetected  had  not  the  remainder  of  the  sky  been  exceptionally  quiet.  Nearly 

the  entire  sky  is  in  medium  or  high  scintillation  on  days  336  and  337,  although 

there  is  a perceptible  tendency  for  lesser  activity  from  -120*  to  -180®  on  336, 

but  from  0®  to  120®  on  337.  This  pattern  is  consistent  with  the  maximum  of  the 

density  structure  passing  over  the  earth  on  day  336,  as  we  noted  from  Figure  11.  , 

In  fact,  comparison  of  Figures  11  and  12  should  make  clear  why  it  is  usually 
necessary  to  use  source  cluster  analysis  when  the  Interplanetary  medium  is  in 
its  more  usually  complex  condition. 

5.  SYNOPTIC  OBSERVATIONS;  SCINTILLATION  SPECTRUM 

The  analysis  of  the  1974  synoptic  index  data  confirmed  our  opinion 
that  it  would  be  highly  desirable  to  have  the  additional  information  available 
which  is  contained  in  the  actual  IPS  spectrum.  We  therefore  analyzed  our  1976 
synoptic  spectral  measurements  for  a signature  of  approaching  turbulent  regions, 
l.e. , a systematic  broadening  of  the  spectra. 

Observations  from  the  University  of  Iowa  COCOA-Cross  radio  telescope 
during  the  summer  of  1976  (17  May-25  August)  were  supplemented  by  observa- 
tions on  the  University  of  Maryland  TPT  radio  telescope,  also  at  Clark  Lake 
Radio  Observatory.  Of  the  151  sources  observed  dally,  a subset  of  41  sources 
proved  favorable  enough  in  location,  flux  and  angular  diameter  to  produce 
measureable  scintillation  power  at  least  once  during  the  observing  period. 

This  subset  of  sources  is  shown  in  geocentric  ecliptic  coordinates  in  Figure 
13.  Of  these, 9 consistently  yielded  spectra  with  sufficient  signal-to-noise 
ratios  to  be  used  in  the  synoptic  analysis  of  spectra  presented  here.  Reasonably 
continuous  coverage  was  obtained  for  the  period  9 June-17  July  (days  166-199)  and 
this  is  the  period  which  will  be  discussed  here. 

The  COCOA-Cross  is  a pencil  beam  instrument  which  operates  at  34.3  MHz 
with  a bandwidth  of  500  kHz.  It  is  electrically  steerable  north-south,  while 
the  rotation  of  the  earth  causes  radio  sources  to  "drift"  through  the  east-west 
antenna  pattern,  providing  usable  data  on  a source  for  about  three  minutes. 

Starting  on  day  183  a rudimentary  east-west  steering  capability  was  implemented 

which  allowed  the  antenna  to  be  pointed  to  three  discrete  east-west  pointings,  i 

1 

symmetrically  spaced  about  the  central  meridian,  and  separated  by  approximately 
two  beam  widths.  This  east-west  steering  was  used  for  sources  3C46,  3C196,  3C216, 
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3C2S4,  3C273,  3C298  and  3C380,  resulting  in  six  to  nine  minutes  of  data  each 
sidereal  day  for  these  sources.  The  latter  four  sources  were  observed  at  two 
of  the  three  possible  east-west  beam  positions. 

The  TFT  was  operated  In  the  pencil  beam  mode  at  a frequency  of  38.0 
MHz  and  with  a bandwidth  of  either  150  kHz  or  800  kHz  depending  on  the  zenith 
angle  of  the  source.  Since  the  TFT  Is  fully  steerable  In  all  directions  (but 
In  discrete  Increments) , It  was  possible  to  track  a source  for  an  extended 
period  of  time  centered  on  central  meridian  transit.  A typical  observation 
consisted  of  Incrementally  tracking  a source  (similar  to  the  east-west  track- 
ing by  the  COCOA-Cross  described  above)  for  up  to  one  hour.  After  allowance 
for  the  beam  shape  this  results  In  approximately  25  minutes  of  data  which  can 
be  used  for  spectral  analysis. 

The  analog  signals  from  both  radio  telescopes  were  passed  through  4 
pole  low-pass  filters  (which  were  designed  to  eliminate  the  effects  of  aliasing) , 
digitized  at  a rate  of  10  samples  per  second,  and  stored  on  a digital  tape 
recorder.  The  observations  from  the  two  radiotelescopes  frequently  overlapped 
In  time  and  generally  showed  good  agreement  between  their  respective  power 
spectra.  Therefore,  the  small  difference  In  observing  frequency  between  the 
two  arrays  was  Ignored  In  the  analysis  and  the  two  data  sets  merged  Into  one 
coaq>o8lte  set  consisting  of  the  highest  quality  observations  frcxn  each  array. 


Data  Reduction 

The  power  spectrum  analysis  consists  of  first  taking  out  long-term 
trends  by  passing  the  data  through  a digital  high-pass  filter  with  a cutoff 
frequency  of  0.02  Hz.  These  high-pass  data  are  then  broken  Into  256  point 
segments  (75.6  sec)  and  a fast  Fourier  transform  Is  performed  on  each  segment. 
The  average  magnitude  squared  of  each  Fourier  component  Is  then  calculated  and 
corresponding  components  from  each  256  point  transform  are  added  together  to 
give  the  desired  spectral  estimates.  On  those  sources  which  are  observed  more 
than  once  per  day  all  of  the  spectral  estimates  are  further  combined  Into  one 
composite  spectral  estimate  with  Improved  stability.  In  order  to  further 
Increase  the  statistical  stability  of  the  spectral  estimates,  the  spectra  are 
smoothed  with  a Hanning  spectral  window.  For  example,  for  the  spectrum  shown 
In  Figure  14a,  which  represents  9 minutes  of  data,  there  are  112  degrees  of 
freedom. 


Two  curves  are  plotted  In  Figure  14a.  The  upper  curve  Is  the  raw 
uncorrected  spectrum  which  consists  of  a scintillation  component  plus  a white 
noise  component  due  to  system  noise.  The  actual  data  points  are  at  the  left 
hand  end  of  each  horizontal  line.  The  sharp  cut-off  above  3 Hz  Is  due  to  the 
4 pole  low-pass  filter  mentioned  above.  Since  In  our  observations  the  scin- 
tillation component  of  the  spectrum  has  always  been  found  to  be  negligible 
above  2 Hz,  we  estimate  the  white  noise  amplitude  by  averaging  across  the  2 
to  3 Hz  range  and  subtract  that  estimate  from  each  of  the  spectral  points. 

This  yields  the  lower  curve  which  for  convenience  has  been  renormalized. 

When  displayed  on  a log-log  plot,  the  spectra  typically  displayed  a 
flat  plateau  from  very  low  frequencies  turning  over  Into  a power-law  at  high 
frequencies,  the  turn-over  occurring  at  anywhere  from  less  than  0.1  Hz  up  to 
0.5  Hz.  Both  the  turn-over  frequency  and  the  rate  at  which  the  spectrum  turns 
over  Into  the  powr-law  provide  Information  about  the  velocity  of  the  scatter- 
ing medium  and  Its  density  turbulence  profile.  The  slope  of  the  power-law  also 
provides  information  about  the  turbulence  profile.  Therefore,  the  spectra  are 
parametrized  by  recording  the  power-law  slope  (a) , the  frequency  at  which  the 
smoothed  spectrum  drops  to  3 dB  below  the  plateau  (v^) , and  the  frequency  at 
which  It  drops  to  6 dB  below  the  plateau  (v^) . These  three  parameters  then 
serve  as  useful  characterizations  of  the  spectra  for  the  analysis  of  large 
amounts  of  data.  We  are  principally  concerned  with  and  In  this  report; 
or  will  not  be  discussed  here,  but  was  analyzed  In  the  paper  by  Mitchell  (1977). 

We  have  taken  several  steps  to  cast  out  observations  that  are  conta- 

minated by  Interference  of  terrestrial  origin.  The  10  sample  per  second  high- 
pass  filtered  data  from  which  the  power  spectrum  was  calculated  is  plotted  for 
every  observation.  Having  Inspected  similar  plots  for  a number  of  cases  of 
known  Interference,  It  was  frequently  possible  to  Identify  terrestrial  Inter- 
ference In  the  observations.  To  further  validate  the  measurements  we  plot  2.5 

sec  and  50  sec  averages  of  the  main  beam  response  and  the  fluctuation  power  in 

the  0.1  and  1.5  Hz  range  (Figure  14b).  The  peak  In  scintillation  power  must 
approximately  coincide  with  the  peak  In  the  main  beam  response,  and  on  the  fine 
time  scale  the  fluctuations  In  scintillation  power  must  be  reasonably  distributed 
throughout  the  passage  of  the  source  through  the  beam,  or  the  observation  will 
be  rejected.  Finally,  the  spectra  themselves  occasionally  show  strong  harmoni- 
cally related  components  due  to  interference,  although  usually  In  this  case  one 
of  the  other  tests  also  confirms  Its  presence. 
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with  these  strict  tests»  we  estimate  that  approximately  two  percent 
of  the  remaining  observations  are  unknowingly  contaminated  by  Interference. 
However,  because  of  the  large  number  of  observations  analyzed,  the  contamina- 
ted observations  should  have  a negligible  effect  on  the  final  conclusions. 

Data  Analysis 

It  has  been  shown  by  Mitchell  (1977)  that  IPS  observations  taken  at 
34  MHz  at  large  elongation  cannot  be  modelled  by  a single  thin  phase  changing 
screen  on  which  most  scintillation  work  Is  based.  The  question  then  arises  as 
to  what  changes  might  occur  In  the  shape  of  the  spectrum  as  the  dominant  scat- 
tering region  (due  to  the  enhanced  turbulence  associated  with  a solar  wind 
stream),  approaches  the  earth.  Having  analyzed  hundreds  of  spectra  we  can  say 
that  the  form  of  the  spectrum  (not  Its  height  or  breadth)  Is  approximately 
similar  when  plotted  on  a log-log  plot.  However,  the  spectra  tend  to  turn 
over  more  sharply  as  they  get  broader.  We  Interpret  this  trend  to  be  a result 
of  the  thickness  of  the  scattering  medium.  In  the  weak  scattering  Interpreta- 
tatlon,  the  superposition  of  contributions  along  the  llne-of-slght  deforms  the 
spectral  shape  of  the  spectrum  from  that  of  a thln-screen  spectrum  and  produces 
a final  spectrum  with  a shallower  slope  and  more  gradual  turn-over  than  any  of 
the  contributing  thln-screen  spectra  would  have.  For  broad  spectra,  the  rela- 
tively near-earth  medium  dominates  the  scattering,  producing  an  approximately 
thin  screen  spectrum.  The  narrow  spectra  are  apt  to  be  the  result  of  scatter- 
ing from  an  extended  region  (for  elongations  > 60°)  In  which  contributions  from 
the  distant  medium  add  low  frequency  power  to  the  spectrum. 

There  are  occasional  spectra  which  differ  significantly  from  the 
majority.  The  most  common  deviant  spectra  exhibit  an  enhancement  at  low  fre- 
quencies (<(  0.2  Hz),  flatten  Into  a well  defined  plateau  at  Intermediate  fre- 
quencies, and  turn-over  again  at  higher  frequencies. 

Dally  values  of  ^or  a number  of  radio  sources  for  the  period  of 
days  166-199,  1976  are  shown  In  Figure  15.  We  have  not  Included  the  two  compo- 
nent spectra  In  this  figure,  since  they  occurred  on  only  a small  number  of  days. 
For  comparison  we  also  display  the  solar  wind  velocity,  density,  and  most  proable 
thermal  speed  measured  by  the  MIT  plasma  instrument  on  the  IMP-7  and  8 spacecraft 
(Solar  Geophysical  Data).  A number  of  cases  of  spectral  broadening  can  be  seen, 
some  of  which  are  correlated  with  solar  wind  streams  as  revealed  by  spacecraft 
data.  There  Is  no  reason  to  believe  that  any  of  these  broadenlngs  are  due  to 
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flare  associated  blast  waves,  since  solar  activity  was  at  a very  low  level 
during  this  entire  period.  In  fact,  as  reported  In  Solar  Geophyaioal  Data, 
for  both  June  and  July  there  were  no  observed  Type  II  or  Type  IV  radio  bursts, 
no  confirmed  H^flares  greater  than  sub  flare  Importance,  and  no  high  energy 
electron  or  proton  events. 

The  spectral  measurements  In  Figure  15  show  a definite  broadening 
both  before  and  during  the  rise  of  the  high  speed  stream  of  day  182.  Since 
our  data  coverage  was  best  for  this  event  we  will  discuss  It  first  and  then 
conment  on  other  prominent  features  for  which,  unfortunately,  we  had  less 
continuous  coverage. 

Day  182  Solar  Wind  Stream 

A period  of  enhanced  IPS  activity  began  late  on  day  181.  This 
occurred  several  hours  before  the  solar  wind  velocity  measured  at  earth  started 
to  rise  due  to  the  approaching  corotating  high  speed  stream,  and  at  a time  when 
the  solar  wind  density  at  earth  had  already  reached  an  enhanced  level.  In  order 
to  study  the  temporal  and  spatial  evolution  of  the  scintillation  power  spectrum 
we  plot  In  Figure  16  quantized  values  of  In  geocentric  ecliptic  coordinates 
with  the  sun  located  at  the  top  of  each  panel.  The  sequence  of  observations 
begins  at  the  sun  and  runs  In  a counterclockwise  direction.  Figure  16a  shows 
that  all  sources  for  which  we  had  data  displayed  spectra  of  average  width.  In 
Figure  16b,  the  first  source  observed,  3C196,  continues  this  trend  and  shows 
no  signs  of  detecting  the  approaching  corotating  solar  wind  stream.  Approxi- 
mately one  half  hour  later  the  situation  changes  significantly  when  3C216  shows 
a broad  spectrum  which  reflects  the  approach  of  the  oncoming  stream.  The  dif- 
ferent response  of  3C196  and  3C216  Implies  a large  spatial  gradient  In  the 
strength  of  the  turbulence  of  this  stream  out  of  the  plane  of  the  ecliptic. 

All  other  sources  shown  In  Figure  16b  displayed  unusually  broad  spectra  except 
3C263.1  which  had  a two  component  spectrum,  the  second  component  of  which  was 
broad.  Figure  16c  shows  the  return  of  scintillation  activity  to  normal;  while 
3C196  Is  broadened,  3C216  and  3C254  have  returned  to  their  average  width.  Con- 
fusing the  simple  picture  of  the  event  terminating  as  abruptly  as  It  had  begun 
is  the  enhanced  width  of  the  eastern  radio  source  3C270.1.  This  behavior  may 
simply  Indicate  a gralnlness  In  the  turbulence  out  of  the  plane  of  the  ecliptic. 
In  any  case  it  Is  clear  from  Figure  IS  that  the  scintillation  spectra  of  all 
other  sources  return  to  average  or  below  average  width  approximately  by  the  end 
of  day  183. 
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Another  aspect  of  the  scintillation  activity  on  day  182  which  Is 
not  brought  out  In  Figure  16  is  the  large  number  of  usually  non-sclntlllatlng 
sources  (at  34  MHz)  that  scintillated  on  that  day.  Because  they  were  rarely 
observed  to  scintillate  we  were  unable  to  classify  their  spectra  based  on 
their  own  past  performances.  This  "all-sky"  enhancement  of  IPS  In  associa- 
tion with  an  enhancement  In  solar  wind  density  (which  precedes  the  velocity 
maximum)  makes  It  clear  that  the  event  under  discussion  Is  similar  to  the 
"all-sky"  events  discussed  by  Erskine  et  at.  (1977). 

Day  169  Solar  Wind  Stream 

Solar  wind  proton  density  peaked  at  approximately  0600  hr  UT  on  day 
169.  There  were  no  radio  observations  on  day  168.  The  only  observation  during 
a 24-hour  period  centered  on  169/0600  is  that  of  4C+21.53  at  169/0940.  This 
source  did  not  show  the  expected  broadening  but  due  to  Its  high  ecliptic  lati- 
tude (42*)  and  large  elongation  angle  at  this  time  (127*)  it  would  be  unwar- 
ranted to  conclude  that  this  one  observation  contradicts  the  enhanced  spectral 
broadening  associated  with  the  solar  wind  stream  described  above. 

Broadening  of  3C48  Spectrum  on  Day  191 

On  day  191  the  scintillation  spectrum  of  3C48  broadened  almost  as 
much  as  It  had  during  the  "all-sky"  event  of  day  182.  The  solar  wind  data 
available  to  us  had  a large  data  gap  which  was  not  over  until  the  middle  of 
day  191  (at  approximately  the  time  when  3C48  was  observed) . Based  on  this 
limited  data  set  there  Is  no  evidence  of  a high  speed  solar  wind  stream  en- 
gulfing the  earth.  Any  enhanced  turbulence  must  have  been  confined  to  the 
northern  ecliptic  hemisphere  at  a distance  no  greater  than  several  tenths  of 
an  AU  above  the  earth.  This  out  of  the  ecliptic  turbulence  gave  rise  to 
broadened  spectra  for  3C48  and  3C196,  and  two  component  spectra  for  3C68.2 
and  3C254.  It  is  unfortunate  that  our  near-ecllptlc  radio  source  3C273  was 
not  successfully  observed  due  to  terrestrial  Interference.  Judging  from  the 
lack  of  a solar  wind  stream  at  earth  we  presume  that  the  spectrum  of  3C273 
would  have  been  of  average  width. 

Day  197  Solar  Wind  Stream 

Radio  frequency  Interference  severely  limited  our  coverage  of  the 
scintillation  activity  associated  with  this  solar  wind  stream.  Since  the 
well  defined  proton  density  enhancement  is  approximately  centered  on  the 
middle  of  day  197t  we  might  expect  Che  associated  scintillation  activity 
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to  be  wholly  confined  to  day  197.  Sources  3C273  and  3C298  both  exhibit 
broadened  spectra  but  3C196  and  3C270. 1 had  spectra  of  average  width.  The 
proton  density  at  earth  had  already  returned  to  normal  when  3C196  (which  was 
in  the  east)  was  observed  at  20  hr  UT  on  day  197.  Its  average  response  Is 
therefore  understandable  since  the  llne-of-slght  did  not  Intersect  the  denser, 
presumably  more  turbulent,  material.  Thus  two  out  of  the  three  sources  which 
were  expected  to  broaden  actually  did  so.  The  one  source  which  did  not  (3C270.1) 
may  have  failed  to  do  so  because  of  a latitudinal  structure  In  the  solar  wind 
turbulence  since  3C270.1  Is  at  a higher  ecliptic  latitude  than  either  3C273  or 
3C298  (see  Figure  13) . 

It  Is  also  Interesting  that  3C273  and  3C298  were  somewhat  broadened 
on  day  196.  Since  these  sources  were  favorably  located  In  the  east  It  appears 
that  they  were  responding  to  the  approaching  corotating  stream  more  than  one 
day  before  It  reached  earth.  The  spectra  for  3C273  on  days  195,  196  and  197 
are  shown  In  Figure  17. 

Interpretation 

We  have  established  In  the  preceding  section  that  the  IPS  spectrum 
broadens  as  a corotating  solar  wind  stream  sweeps  over  or  very  near  to  the 
earth.  We  can  discuss  this  phenomenon  qualitatively  In  terms  of  the  usual 
thin  phase  changing  screen  model  as  the  source  of  the  scintillations.  Mitchell 
(1977)  has  shotm  that  IPS  observations  at  34  MHz  require  a thick  screen  model, 
but  If  the  effects  are  similar  for  the  screens  whose  contributions  dominate 
the  spectrum,  the  thin  screen  model  should  give  adequate  insight  Into  the 
spectral  behavior.  The  high  frequency  dependence  of  our  best  spectra  for 
3C48  can  be  characterized  by  a power-law  so  we  shall  Interpret  our  results 
assuming  that  the  spatial  turbulence  spectrun  was  also  a power  law. 

There  are  two  possible  mechanisms  which  control  the  spectral  width 
from  a given  thin  screen  which  can  be  discussed  In  terms  of  the  theoretical 
concepts  presented  earlier  In  Section  1.  In  the  first,  for  a point  source, 
or  for  a very  small  apparent  diameter  source  such  as  3C48,  Fresnel  diffraction 
"filters  out"  much  of  the  energy  associated  with  the  large  scale  turbulence 
and  causes  the  temporal  power  spectrum  to  turn  over  at  V^/^^^where 
Is  the  bulk  velocity  of  the  turbulence  perpendicular  to  the  llne-of-slght, 
and  Z Is  the  distance  from  the  observer  to  the  scattering  region.  The  second 
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mechanism  (finite  angular  size  cutoff)  applies  to  larger  diameter  sources 
such  as  3C273,  where  the  turn-over  frequency  Is  Vc  V /Z.  Thus  regardless 
of  which  mechanism  Is  dominant.  Increasing  the  perpendicular  component  of 
the  bulk  velocity  or  decreasing  the  distance  to  the  scattering  region  causes 
the  spectrum  to  broaden. 

In  order  to  settle  the  question  of  whether  distance,  velocity,  or 
both  were  the  cause  of  the  observed  broadening  we  examined  the  solar  wind 
velocity  and  spectral  data  obtained  by  multi-size  IPS  observations  at  74  MHz 
by  the  UCSD  group  (<7.  Harmon^  Pipivate  Camunioation) . Due  to  randomly  spaced 
gaps  In  our  observations  and  limited  overlap  in  the  available  data  sets,  we 
were  left  with  only  a sequence  of  observations  of  3C48  centered  on  day  182 
which  could  be  compared  at  both  frequencies.  The  velocity  derived  from  the 
74  MHz  observations  of  3C48  showed  no  Increase  on  day  182  or  for  two  days 
thereafter.  On  the  other  hand  both  the  34  and  74  MHz  IPS  spectra  were  anom- 
alously broadened,  leading  us  to  the  conclusion  that  at  least  for  this  one 
well  documented  case  and  for  the  assumption  of  a power-law  turbulence  spec- 
trum It  Is  the  changing  distance  to  the  turbulent  region  which  controls  the 
width  of  the  spectrum.  Considering  the  all-sky  nature  of  the  day  182  IPS 
activity  this  Is  In  agreement  with  the  conclusion  of  Erakine  et  aZ.  (1977) 
that  all-sky  events  are  not  a response  solely  to  Increased  velocity,  but  are 
associated  with  enhanced  turbulence  near  the  leading  edge  of  a high  speed 
stream  where  the  density  Is  also  enhanced. 

6.  CONCLUSIONS 

Our  analysis  of  the  May-December  1974  synoptic  IPS  Index  data  has 
led  ts  to  the  conclusion  that  source  sizes,  typically  ^ one  arc-second  at  34.3 
MHz, often  limit  the  IPS  response  from  those  sources  to  within  a few  tenths  of 
an  AU  from  earth  If  there  Is  enhanced  turbulence  with  the  vicinity  of  earth. 
Comparison  with  spacecraft  solar  wind  density  measurements  shows  that  local 
high  densities  lead  to  strong  IPS  response.  On  the  other  hand,  if  densities 
near  earth  are  low,  there  may  be  IPS  response  to  more  distant  turbulence, 
especially  for  elongation  angles  c < 90* . 

Our  May-August  1976  synoptic  IPS  spectrum  data  Indicate  that  the 
width  of  the  scintillation  power  spectrum  can  provide  Information  on  the  diatccnoe 
of  the  turbulent  region  from  the  earth,  within  the  range  which  can  be  probed  due 
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Co  Che  flnlCe  source  specCral  filtering  effect.  Progressive  broadening  of  specCra 
were  observed  as  corocaclng  curbulence  approached  Che  earch. 

In  nelcher  1974  or  1976  were  chere  sufflclenCly  sCrong  flare-assoc laced 
solar  wind  shocks  to  determine  Che  IPS  response  at  34.3  MHz  to  these  structures 
which  are  more  common  during  the  maximum  of  the  solar  cycle.  However,  the  pro- 
nounced recurrent  solar  wind  structure  during  the  decline  of  the  last  cycle 
allowed  us  to  establish  Chat  the  IPS  reponse  signature  was  clearly  associated 
with  the  enhanced  densities  in  the  leading  edges  of  the  corotaClng  screams  which 
should  be  similar  to  Chose  In  an  expanding  flare-associated  shock. 

A nominal  range  of  several  tenths  of  an  AU  for  probing  the  Interplan- 
etary medium  using  UPS  observed  at  34.3  MHz  means  that  reliable  prediction  may  be 
feasible  If  a lead  time  of  about  a day  before  the  onset  of  a geomagnetic  storm 
is  acceptable.  This  time  Interval  should  be  about  the  same  for  either  corota- 
ting streams  or  flare-produced  shocks.  Of  course.  It  is  necessary  to  have 
sources  near  the  ecliptic  plane  for  positive  Identification  of  turbulence 
about  to  envelope  the  earth.  One  must  bear  in  mind  that  Che  diatrlbutlon  of 
scintillating  sources  is  not  unifoirm  over  the  sky  (see  Figure  13).  The  orienta- 
tion of  the  galactic  plane  Is  such  that  as  seen  from  California  the  source  grid 
is  densest  toward  the  east  in  Spring  and  Summer  and  sparsest  In  the  Fall.  The 
sensitivity  to  flare  shock  disturbances  will  likewise  be  least  In  early  Winter 
therefore  prediction  reliability  will  depend  somewhat  on  time  of  year,  although 
the  least  favorable  periods  for  either  corotating  or  flare-associated  distur- 
bances last  less  than  two  months. 

One-day  prediction  capability  can  only  be  established  with  high  quality, 
continuous  synoptic  data.  We  were  unable  to  obtain  data  meeting  these  specifi- 
cations during  the  period  of  this  contract.  Therefore  our  results  are  not  definitive 
however  we  believe  that  we  can  draw  the  following  conclusions: 

Prediction  of  the  onset  of  geomagnetic  disturbances  using  IPS  observa- 
tions at  34.3  MHz  does  not  appear  to  be  reliable  for  lead  times  exceeding  one 
dayt  but  there  is  preliminary  evidence  that  it  may  be  feasible  for  < one-day  lead 
times,  especially  if  synoptic  spectral  observations  are  utilized. 
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Thln-scrMO  •paetruu  calculated  for  five  distances  L (out  along 
the  line  of  sight)  by, using  expansions  Jj*  Jix>  A&d  Jxii«  Their 
relative  nonalisatlon  Is  given  by  F « r~^i  although  the  absolute 
noraallsatlon  of  F^  Is  arbitrary.  The  Integrated  eclntUlatlon 
power,  integrated  along  the  line  of  sight  for  each  frequency  with 
F^  « r~  , Is  given  by  the  heavy  solid  line,  (a)  A relatively 
"swell**  source  radius  Q ■ 0.1**.  (b)  A larger  source  ■ 0.5", 

dewonstratlng  the  strong  suppression  of  the  high-frequency  power 
fwoa  wore  distant  screens. 
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Figure  2 Contribution  to  the  scintillation  Index  (band**passed  from  0.1  to 

1.5  Hz)  as  a function  of  distance  along  the  line  of  sight  for 

three  source  angular  diameters  (6).  The  calculation  vas  done  for 
_4 

F « r • q ~ 3.6,  solar  wind  velocity  V ■■  400  km/sec,  and  34.3  MHz. 


Figure  3 
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GeoMtry  of  a corotating  turbulent  region  with  varloua  relative 
poaltlona  of  the  earth  and  the  llnea  of  al^t  to  aourcea  at  60* 
and  120*  elongation.  The  locatlona  of  the  centers  of  slabs  7 
and  13  (of  the  total  of  20)  are  Indicated. 


1974  DOY 

Figure  4 Synoptic  observations  of  6 scintillating  sources  In  s region  to 
the  west  of  the  sun  of  dlasnslons  f*  (ecliptic  latitude)  by  57* 
(ecliptic  longitude).  Also  shown  ere  the  following:  IMP-7  end 
8 solar  wind  proton  density  velocity  (Los  Alaaos  Scientific 
Laboratory);  Gsonsgnstlc  A^  Index,  and  the  occurrences  of  gso- 
■sgnstlc  sudden  coasenesaants  (SSC)  and  Type  II  and  IV  solar 
bursts  In  close  association  [Solar  Gsophyslcsl  Data];  and  Inferred 
Intsrplanstary  asgnstlc  field  direction  [Svalgaard,  1975]. 
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Flgur*  5 l«dlo  MurcM  displaying  uMbla  IPS  Indax  during  Hny-Daeaad>ar  1974 
elustarad  for  IPS  sky-box  analyala.  For  this  gaocsntrlc  all-sky 
pro j action  tba  aslauthal  coordlnata  Is  acllptle  longltuda  ralatlva 
to  tha  Vamal  Equinox  (dagraas)  and  ths  radial  coordlnata  Is  acllptle 
latitude  (dagraaa). 
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Figure  6 Rasulce  of  radio  aourca  clustarlng  analyals  for  1974,  daya  150-190. 

Shading  la  aach  aky-boz  Indlcataa  high  (dark  ahadlng),  aodarata 
(aadlia  ahadlng),  low  (light  shading),  or  alaalag  (no  shading)  IPS 
Indaz.  Haavy  llnas  running  diagonally  In  tha  body  of  tha  plot  of 
acllptlc  longitude  va.  tlae  rapraaant  tha  position  of  tha  sun  and 
tha  dlraetlon  180*  froa  tha  sun  (aatl-sun).  Dashed  llnas  running 
diagonally  la  body  of  aach  plot  are  drawn  at  90*  acllptlc  longltuda 
east  and  wast  of  tha  sun.  At  tha  canter  are  plots  of  IPS  'activity 

t***  following  cases t all  10  bozas  averaged  together,  < Q > 
aaatam  sky  boxes  only,  < Q wastam  sky  bozas  only,  < Q > 

high  latltuds  box  9 only,  < Q > . At  tha  bottom  of  tha  ^ot 

tk*  following  corralatlva  datat  IMP-7  and  8 spacecraft  solar  wind 
proton  density  and  velocity;  Interplanetary  aagnatlc  field  esgnltude 
^*zyl  oolar  ecliptic  direction  9^;  Geomagnetic  A^^  Index;  Inferred 
Interplanetary  aagnatlc  field  [Svalgaard,  19751  (dots  with  dashed  llnei 
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rifur*  9 TIm  Mtl**  crocs  eorrslstion  results  for  corrslstlon  of  sll-sky  IPS 
sctlTlty  Indus'  < Q > versus  the  followlnf'  < Q > (eutocorrelstlon)/ 
soler  wind  proton  density  ll»  Interplsnetery  nsgnstlc  field  nsgnltude 
|B  I (ecliptic  projection)*  soler  wind  proton  bulk  velocity  V*  plsne- 
tery  iSfiMUistl  r A Indes*  end  occurrences  of  Interplsnetery  sector 
crossings.  Plotted  Is  tbs  crosscorrslatlon  coefficient  r(T)  for  tine 
shifts  > 10  dey  A r A 159  deys.  The  dete  bese  Is  1974*  days  150-362. 
Also  listed  Is  r (t)  for  the  peek  nearest  t * 0 days*  as  well  as  ths 
nosher  of  points*  n*  In  the  tlae  series.  Ths  diverging  lines  Indicate 
the  99S  confidence  Halts. 
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IPS  *sky  bos'  analysis  and  solar  wind  danaltyN  for  'danalty  avsnt 
of  1974,  day  335. 
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Figure  12 


Equel-eraa  eky  maps  of  IPS  activity  for  1974  a)  day  332 » b)  day 
c)  day  337. 
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Sun,  day  165 


Figure  13 


Subset  of  41  sources  observed  by  the  COCOA-Cross  radio  telsseope. 

Sources  plotted  with  s filled  circle  consistently  yielded  useful  spectra. 
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varlablM  ar*  2.5  Mcond  averages,  and  barred  variables  are  50  second 
triangle  valgbtad  svaragaa. 


I Abrupt  incrsM  in  Kp  # Anomalous  spectrum 


Day 


Figure  15  Top:  veraus  tine  for  9 aourceat  Bottom:  aolar  wind  probe  data 

(Solar  Goophyaioal  DatCt  197?).  The  atart  of  the  three  known  solar 
wind  atreaae  la  highlighted  by  ahedlng.  On  thoae  days  when  a two 
component  spectrum  was  present,  or  the  spectrum  was  ambiguous,  the 
date  la  omitted  end  e diamond  la  plotted  near  the  bottom  of  the 
relevant  curve. 
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Figure  16  (a-c)  for  a number  of  aourcaa,  plotted  in  geocentric  ecliptic 

coordlnatea.  The  "day  cut"  la  at  the  aun,  and  the  obaervatlona  proceed 
In  a countercloclnrlae  aequence.  The  threahold  for  "low"  and  "high"  la  one 
atanderd  deviation  from  tha  average,  (d)  Llne-of-alght  to  3C48  pro- 
jected onto  the  ecliptic  plane  on  1976t  day  182.  Alao  ahown  la  the 
approximate  ahapa  of  tha  canter  of  a corotatlag  denalty  enhancement 
for  daya  181,  182  and  183  ar  tha  time  when  3C48  waa  obearvad. 


ritur*  17  (a)  Saoothad  sclntlllatloii  apactra  for  3C273  on  daya  19S»  196  and  197 

lUuatratlng  prograaalva  broadanlng  aa  a corotating  region  approachaa 
tha  aarth.  (b)  Lina  of  alght  to  3C273  projactad  onto  tha  acllptlc  plana 
on  1976*  day  197.  Alao  ahown  la  tha  approslaata  ahapa  of  tha  cantar  of 
a eorotatlng  danalty  anhancaaant  for  daya  195 » 196  and  197  at  tha  tlaa 
ahan  3C273  «aa  obaarrad. 
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